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ABSTRACT 


Significant advantages of the Variable Polarity Plasma Arc 
(VPPA) welding process include faster welding, fewer repairs, less 
joint preparation, reduced weldment distortion, and absence of 
porosity. The power distribution has been analyzed for an argon 
plasma gas flow constituting the fluid in the VPPA welding process. 
The major heat loss at the torch nozzle is convective heat transfer; in 
the space between the outlet of the nozzle and the workpiece, 
radiative heat transfer; and in the keyhole in the workpiece, 
convective heat transfer. The power absorbed at the workpiece 
produces the molten puddle that solidifies into the weld bead. Crown 
and root widths, and crown and root heights of the weld bead are 
predicted. The basis is provided for an algorithm for automatic 
control of VPPA welding machine parameters to obtain desired weld bead 
dimensions . 
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NOMENCLATURE 


A c = Crown reinforcement area (cm 2 ) 

A e = Cross-sectional area of electrode (cm 2 ) 

A„ = Nozzle circumferential area (cm 2 ) 

A r = Root reinforcement area (cm 2 ) 

C p = Constant pressure specific heat (J/°K) 

D = Nozzle diameter (cm) 

D e = Electrode diameter (cm) 
d c = Crown bead width (cm) 
d r = Root bead width (cm) 
dv 

— = Voltage drop per unit length (V/mm) 
dL 

d v = Feed wire diameter (mm) 
g = Gravitational acceleration (cm/s 2 ) 
h c = Crown height (mm) 

h e xit = Plasma jet enthalpy at the exit of workpiece (J/g) 

h„ = Plasma enthalpy corresponding to the nozzle wall temperature 

(J/g) 

h p u$na = Plasma enthalpy (J/g) 
h r = Root height 

hj # i , h 1/0 = Plasma enthalpy at entrance and exit ends, 

respectively, of gap between electrode and nozzle 

h 2 7 i / h 2 0 = Plasma enthalpy at entrance and exit ends, 
respectively, of nozzle 

h 2 j , h 3 0 = Plasma enthalpy at entrance and exit ends, 
respectively, of standoff column 

h 4 / i, h 4/0 = Plasma enthalpy at entrance and exit ends, 
respectively, of workpiece keyhole 
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I m + , I m “ = Main electric current in straight and reverse 
polarities, respectively 

I + , Ip" = Pilot electric current in straight and reverse 
polarities, respectively 

k = Thermal conductivity of plasma jet 

k E = Thermal conductivity of electrode 

L e = Length of electrode 

= Distance between electrode and nozzle 

L 2 = Axial length of nozzle 

L 3 = Column distance of standoff 

L 4 = Thickness of workpiece 

m = Argon gas mass flow rate 

n = Parabolic exponent of reinforcement shape 

P = External pressure 

P E = Power input at electrode 

Pj j = Joule heating power input in the gap between electrode and 
nozzle 

Pj ' 2 - Joule heating power input in the nozzle 

Pj 3 = Joule heating power input in the standoff column 

P |(1 = Nozzle work function induced power input 

Pjj 2 = Nozzle anode and cathode drops induced power input 

P 2 = Workpiece work function induced power input 

P U/2 = Workpiece anode and cathode drops induced power input 

p i i r p i o = Plasma arc power at entrance and exit ends 
' ' ' ' respectively, of gap between electrode and nozzle 

Pj i , P 2 0 = Plasma arc power at entrance and exit ends 
* ^ respectively, of nozzle 

p 3 i / p 3 o = Plasma arc power at entrance and exit ends 
1 respectively, of standoff column 
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P 4 , i / P 4 t o = Plasma power at entrance and exit ends, 
respectively, of workpiece keyhole 

Q d i s = Plasma arc discharged power loss (at exit end of keyhole) 

Q g = Power loss in the gap between electrode and nozzle 

Q n = Total power loss in the nozzle 

Q„ yC = Convective heat transfer power loss in the nozzle 

Q„ t e = Free electron power loss in the nozzle 

Q„ G = Direct power loss from the surface of the nozzle 

contributed by total power input to work function and part 
of power input to anode and cathode drops at nozzle 

Q s = Standoff column total power loss 

Q v = Total workpiece power loss 

Q H;C = Convective heat transfer workpiece power loss 
Q U/6 = Free electron power loss in the workpiece 
R = Nozzle radius 

Rj = Distance between heat source and point of computation at 
workpiece 

t + , t_ = Time duration for straight and reverse polarities, 
respectively 

T c = Temperature at the point of computation 

T 0 = Environment temperature; Initial temperature of workpiece 
u = Welding traveling speed 

V A , v b / V c , V D , V E , V F , V G = Straight polarity electric voltage, 

shown in Figure 1, at points A, B, 

C, D, E, F, and G, respectively 

V A 7 , V B 7 , V c / , V D 7 , V E 7 , V F 7 , V G 7 = Reverse polarity electric 

voltage, shown in Figure 1, 
at points A 7 , B 7 , C 7 , D 7 , E 7 , 
F 7 , and G 7 , respectively 

• 

V G = Argon plasma volumetric flow rate in the gap between 
electrode and nozzle 

V s = Argon plasma volume in the standoff column 
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V w = Wire feeding speed 

W g = Initial gap distance between two workpieces 
y = Coordinate axis 

y c = Axis of coordinate to be used for spatial temperature 
distribution computation 

Z = Coordinate axis 

Z c = Axis of coordinate to be used for spatial temperature 
distribution computation 

Zi 35 = Heat source spatial point of distribution 
a = Angle; workpiece thermal diffusivity 
B = Angle 

&V a , aV c = Anode and cathode electric potential drops, 
respectively 

AV m + , aV m ~ = Main electric voltage drops in straight and reverse 
polarities, respectively 

aV p + , aV p " = Pilot electric voltage drops in straight and reverse 
polarities, respectively 

c R = Reverse polarity voltage rise 

a = Surface tension* coefficient 

<J) = Fraction of heat source distributed in spatial coordinate 

= Electrode work function 
= Nozzle work function 
0 V = Workpiece work function 

^ = Fraction of nozzle, or workpiece, power loss contributed by 

anode and cathode drops 
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I. 


INTRODUCTION 


The plasma arc is a concentrated energy source commonly used in 
welding and cutting processes. It is composed of a partially ionized 
gas stream produced by forcing an inert gas to flow through an electric 
arc and emerge from a constricting nozzle. With its high energy 
density and velocity, the plasma arc, when impinging on the workpiece, 
can create a hole in the molten pool generated by its heat and can 
penetrate through the workpiece. Depending on the operating 
parameters employed, this hole may either become self-healing or 
remain open as the arc transverses along the workpiece. A "keyhole" 
welding process occurs in the self-healing case, in which the molten 
metal in front of the arc flows around the arc column and resolidifies 
behind the arc. On the other hand, a cutting operation is achieved if 
the hole remains open (O'Brien, 1968) . 

The plasma welding process is essentially an extension of gas 
tungsten-arc welding (GTAW) ; however, constraining the arc to flow 
through a nozzle produces much higher energy density in the arc and 
much higher gas velocity and momentum. The process produces deep 
penetration welds by forming a keyhole in the workpiece by the 
pressure of the gas flow. The metal melted in front of the advancing 
keyhole flows around to the rear where it solidifies to form the weld 
bead. This is distinct from the essentially surface melting produced 
by the GTAW process (which cannot normally penetrate to a depth equal 
to the width of a weld pool) . 

In 1955 Linde Air Products introduced a plasma arc torch for 


6 



metal cutting applications, and by 1965 Linde had developed an 
automatic plasma arc welding facility for Westinghouse Electric 
Corporation to fabricate a 12 0 in. diameter, 3/8 in. thick D6AC steel 
rocket motor case for the Titan III-C booster assembly (Miller and 
Filipski, 1966; Privoznik and Miller, 1966). The plasma arc welding 
process was reported to halve the welding time required. Then in 19 65 
Thermal Dynamics Corporation reported the use of direct current 
reversed polarity (workpiece negative) plasma arc to join 1/4 in. 
thick aluminum plate (Cooper et al . , 1965). At the end of the 1960 ' s , 
Van Cleave at the Boeing Company began his efforts to combine the 
plasma arc process with a variable polarity feature in which the 
electrode polarity was periodically reversed. Alternating 
electrode potential for aluminum welding had been investigated as 
early as 1947 (Herbst, 1947) . Difficulties with welding power 
supplies in this application were evident early-on. In addition, 
when variable polarity plasma arc (VPPA) welding was used in the U.S. 
Army Roland Missile Production Program, development problems such as 
arc pressure pulsation were noted (Nunes et al., 1981; 1983) . As a 
result of Van Cleave's promising work at Boeing, a VPPA research and 
development project was initiated in 1978 at the NASA Marshall Space 
Flight Center, to determine the potential for replacing the GTAW 
system used in the fabrication of the Space Shuttle External Tank 
(ET) . After the original used in the fabrication and process were 
improved, VPPA welding finally exceeded the expectations (Nunes et 
al., 1984a; 1984b). 

The Space Shuttle ET (diameter 28.6 ft., length 154 ft.) is the 
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largest known "drop tank," carrying 140,000 gallons of liquid oxygen 
and 380,000 gallons of liquid hydrogen. From the outset of ET 
production, the conventional GTAW system operating in the direct 
current electrode negative mode was used. With the GTAW system, in 
spite of the requisite careful joint preparation, weld porosity and 
defects had been ever-present problems since the fabrication of the 
Saturn lunar rocket first stage. Therefore, the decision was made at 
the outset to use 100% radiographic inspection. Weld porosity and 
defects had to be ground out and repaired systematically. 

The VP PA process significantly reduced porosity and defects with 
less stringent joint preparation and, in addition, reduced weldment 
distortion, and speeded up the welding process. 

In this study, the power distribution for an argon plasma gas has 
been analyzed along its course through the VPPA welding system. The 
study includes the following sections: 

(1) Electric potential and current distributions of pilot arc 
and main arc circuits during straight and reverse polarity 
processes of VPPA welding system operations; 

(2) Electric power input at electrode, power input in the gap 
between electrode and nozzle, power input within the 
nozzle, power input in the free plasma jet column of 
standoff, and power input in the keyhole of workpiece; 

(3) Thermal power loss in the VPPA welding system including 
electrode heat conduction loss; heat radiation loss within 
the gap between electrode and nozzle; power loss within the 
nozzle which contains heat convection loss within the 
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nozzle, direct power loss from the surface of the nozzle, 
and electron collision power loss within the nozzle; power 
loss in the standoff column; workpiece power loss which 
contains convective heat loss in the workpiece, direct 
power loss from the surface of the keyhole, and electron 
collision power loss in the keyhole; and power discharge 
loss ; 

(4) Plasma arc enthalpy of various sections of VPPA welding 
system including plasma arc enthalpy within the gap between 
electrode and nozzle, plasma arc enthalpy within the 
nozzle, plasma arc enthalpy within the standoff column, 
plasma arc enthalpy within the keyhole of the workpiece, 
and plasma arc enthalpy discharged from the workpiece; 

(5) Prediction of crown and root bead widths of welding; 

(6) Prediction of crown and root heights of welding. 
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II. VARIABLE POLARITY PLASMA ARC WELDING SYSTEM 


ELECTRIC POTENTIAL CALCULATION 


The significant current carriers in a welding arc are electrons 
and positive ions. Electrons carry the bulk of the current, moving 
rapidly from negative cathode to positive anode. The positive ions 
drift more slowly through the interelectrode space. In a fixed 
polarity arrangement, the differential drift rate results in an 
asymmetrical heating of the welding arc's ends. The cathode receives 
less heat and the anode more heat. The straight polarity mode of 
operation entails a negative electrode (cathode) and a positive 
workpiece (anode) . Where the primary object of the weld process is to 
deliver the maximum heat to the workpiece with minimum deterioration 
of the electrode, straight polarity is used. 

Reverse polarity, which entails positive electrode and negative 
workpiece, has the advantage that the workpiece is subjected to a 
cleaning process (cathodic cleaning) by the impingement of positive 
ions on the workpiece surface. Recent work suggests that heavy 
positive ion weight is irrelevant to cleaning process which seems to 
be more like capacitor breakdown than hammer impact. In the case of 
plasma arc welding, reverse polarity action appears to condition the 
surface of the aluminum alloy so that the molten metal flows easily and 
controllable under the arc. It is conjectured that this fluid 
control is accomplished through breaking up of surface oxide films. 
Cutoff of reverse polarity during VPPA welding transforms a weld metal 
flow which closes smoothly and soundly behind the keyhold to an 
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irregular, intermittent, globular flow. This process leaves a 
rough, lumpy bead sunken below the parent metal surface and protruding 
jaggedly from the root of the weld. However, continuous reverse 
polarity is not necessary to provide adequate cathodic cleaning 
action. 

The variable polarity square-wave with unequal straight and 
reverse polarity time intervals offered a combination of the high 
heating capability of straight polarity with the cleaning feature of 
reverse polarity. Adequate cleaning is obtained by incorporating a 
relatively short (1/10 to 1/5 the duration of the straight polarity 
current) pulse of reverse polarity current into the welding current 
waveform. 

Figure 1 show the VPPA electric power system. During the 
straight polarity mode of operation, it employs a negative electrode 
(cathode) and a positive workpiece (anode) ; during the reverse 
polarity mode of operation it employs a positive electrode and a 
negative workpiece. Figure 1 at left illustrates the straight 
polarity mode; at right, reverse polarity mode. 

In electric potential calculation, one has to consider 
electrode, nozzle and workpiece work functions; anode and cathode 
drops; electric potential intensity, etc. In reverse polarity, in 
addition to those considered earlier, one has to consider extra 
electric potential rise to compensate the thermionic emission in 
metal workpiece which is very much similar to the extra cathode drop. 
Following these logics, calculation of electric potential can be 
obtained from the following formulations: 
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(A) For Straight Polarity 

In straight polarity, electron streams flow from plasma arc to 
workpiece (anode) . Summation of workpiece work function and anode 
drop is balanced by the electric potential difference between anode 
potential V F and plasma arc potential V E , i.e., 

V F - V E = AV a + * v 
or 

V E = V F - AV a - * w 

It is grounded at point F which implies that V F = 0. 

Electron traveling distance between points E and C is L 2 + L 3 
which means that electric potential difference between these two 
points along the direction of electron flow is 

dV 

V E - V c = — (Lj + L 3 ) 
dL 

or 

dV 

V c = V E (L 2 + L 3 ) 

dL 

Similarly, electric potential difference between V„ and V c 
resemble to the electric potential difference between V F and V E , shown 
earlier, which is 

v D - v c = AV a + 4> x 

or 

V D = V c + aV a + t. 

Electric potential difference between V c and V B also resembles 
to the electric potential difference between V E and V c , namely, 
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dV 

V C - V B = — L x 
dL 

dV 

V B = V c L x 

dL 

Electric potential difference between cathode electric 
potential V A and the neighborhood plasma arc potential V B , V A - V B is 
balanced by the electrode work function and negative cathode drop, 
i.e. , 

V* - V B = * E - AV C 
or 

v A = v B - AV C + * E 

The total potential drop in main arc aV n + is V F - V A - - V A . Thus, 
4V„* * - V, 

The pilot arc potential drop &V p + is V D - V A , i.e., 

= V D - v A 

where <p z = Electrode Work Function 
= Nozzle Work Function 
= Workpiece Work Function 
AV a = Anode Drop 
aV c = Cathode Drop 
dV 

— = Electric Potential Intensity 

dL 

(B) For Reverse Polarity 

In reverse polarity, similar calculation can be accomplished for 
solving the electric potential differences. To give an example, the 
electric potential difference between cathode electric potential V F ' 
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and the neighborhood plasma arc potential V E ' resembles to the 
electric potential difference between V A and V B in straight polarity, 
in addition to the consideration of extra reverse polarity rise, e R , 
to be used to compensate the thermionic emission loss in workpiece 
metal (#2219 Aluminum), namely 

V E # = V F ' “ + aV C + «R 

In this study, mathematical expression based on Nunes (1991a) has been 
adopted, i.e., 

t_ 

e R = Reverse Polarity Rise =8+20 

t + + t + 

The rest of the mathematical expressions of electric potential 
f erences in reverse polarity can follow the similar logics shown in 
the calculation of straight polarity. These mathematical 
expressions are shown as follows: 

dV 

V c ' = V E ' + — (Lj + L 3 ) 
dL 

dV 

Vp 7 — V c ' — — L 2 + AV a + (J)^ 
dL 

dV 

v b' “ V c ' + — L t 
dL 

v a 7 = V B ' + aV a + # E 

= aV A , y , = V A , - Vp, = Vj , 

* V p" = * V A'D* - V D . - V*. 

where t + = Straight Polarity Time Period 
t_ = Reverse Polarity Time Period 
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III. VPPA Welding System Electric Power Input 


The pilot arc current, superimposed upon the main arc current 
between the electrode and the plasma nozzle, comes from a separate 
power supply located in the plasma control console. The pilot arc 
polarity does not alternate. The nozzle is held positive with 
respect to the electrode so that the electron flow is always away from 
the electrode toward the nozzle. The main arc is started by applying 
high frequency alternating voltage between the electrode and the 
plasma nozzle so as first to establish the pilot arc from which the 
main arc proceeds. The high frequency current is controlled by the 
plasma control console, but is generated by the spark gap oscillator 
located in the hot block. 

The total electric power input of the VPPA welding system 
includes the pilot arc between electrode and nozzle as well as the main 
power arc between electrode and workpiece during the straight 
polarity and reverse polarity modes. VPPA Electric Power Input can 
be calculated as follows: 

Im + *Vt + + Ip + *Vt + + + Ip’AVp-t. 

Ptotal = t 3 " 1 ) 

t + + t_ 

where I m + = Straight Polarity Main Current 
I p + = Straight Polarity Pilot Current 
I m " = Reverse Polarity Main Current 
I p " = Reverse Polarity Pilot Current 
Equation (3-1) is the total electric power input of the system which is 
the summation of power input for the plasma arc jet in the gap between 
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electrode and nozzle, power input within nozzle, power input in free 

jet column (standoff) , and power input within keyhole of workpiece. 

These power input can be itemized, and calculated as follows: 

(III-A) Power Input at the Electrode 

(i m + + Ip + )(AV C - * E )t + + (I m -- i p -)(AV a - * E )t_ 

p £ = (3-2) 

t + + t. 

(III-B) Power Input in the Gap Between Electrode and Nozzle 
(Joule Heating) 


J / i 



f dv 1 


f dV 1 


+ Ip + ) 

— L, 

t + + (I*’ - Ip") 

— Li 

t 


IdL J 


IdL J 



(3-3) 


t + + t. 


(Ill— C) Power Input Within Nozzle 

This includes power input through the work function of metal 

+ v*„t_ 


K / 1 


(3-4) 


t + + t. 


Power input through anode and cathode drops, 
I p + AV a t + + Ip-AV.t. 


« / 2 


(3-5) 


t + + t. 


and power input through Joule heating of nozzle 


(dV 


dLj 


(dV 


L 2 t + + (I,- - Ip-) 


VdL 


L 2 1 _ 


J , 2 


(3-6) 


t + + t. 


(III-D) Power Input in the Free Plasma Jet Column (Standoff) 
(Joule Heating) 
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L 3 t + + I n 


t + + t. 


L 3 t_ 


(3-7) 


(III-E) Power Input in the Keyhole of Workpiece 


This includes power input through the work function of metal and 


reverse polarity potential rise 

^w^-+ + ( £ R ” ^w)^- 

P v,l = 

t + + t_ 

and power input through anode and cathode drops 
l/av a t + + I„-AV c t_ 

Pv , 2 

t + + t _ 

(II - F) Total Power Input of VPPA Welding System 


( 3-8a) 


( 3 -8b ) 


Total power input of VPPA welding system, shown in Equation (3-1) 
is equal to the summation of power input shown in Equation (3-2) to (3- 
8) , namely 


total P E +P J , 1 +P X , 1 +P M , 2 +P J , 2 +P J , 3 +P v , 1 +P v , 


(3-9) 
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IV. Power Loss Calculation of VPPA Welding System 


Power loss of VPPA system includes (a) electrode heat conduction 
loss, (b) heat radiation loss within the gap between the electrode and 
the nozzle, (c) power loss within the nozzle which contains heat 
convection loss within the nozzle, direct power loss from the surface 
o.f the nozzle, and electron collision power loss within the nozzle, 
(d) power loss in the standoff column, (e) workpiece power loss which 
includes convective heat loss in the workpiece, direct power loss from 
the surface of the keyhole, and electron collision power loss in the 
keyhole, and (f) power discharge loss. These power losses of VPPA 
welding system are itemized and are illustrated as follows; 

(IV-A) Electrode Heat Conduction Loss 

One-dimensional formulation of electrode heat conduction loss 
is used. The equation is given by 

(T e - T 0 ) 

Qe = k E • A e (4-1) 

Le 

where 

Q e = Electrode Heat Conduction Loss 
k E = Thermal Conductivity Coefficient of Electrode 
A e = Cross-sectional Area of Electrode 
L e = Length of Electrode 

T e - T 0 = Temperature Difference between Ends of Electrode 
(IV-B) Heat Radiation Loss Within the Gap Between Electrode 
and Nozzle 

Plasma arc heat radiation loss within the gap between electrode 
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and nozzle can be computed from the following equation (Evans and 
Tankin, 1967; Cram, 1985; and Nunes, 1991b): 


Q g (W) = V G (m 3 ) (4 x 10 10 


T(°K) 

15000 



T(°K) 


1 6 


15000 


(4-2) 


where 

Q g = Heat Radiation Loss Within the Gap Between Electrode 
and Nozzle (W) 

V G = Plasma Arc Jet Volume Within Gap Between Electrode 
and Nozzle (m 3 ) 

T = Average Enthalpy of Plasma Arc Jet Corresponding 
Temperature (°K) 

(IV— C) Power Loss Within the Nozzle 

Power Loss within the nozzle is comprised in the following items: 
(IV-C-1) Heat Convection Loss Within the Nozzle 

Heat convection loss within the nozzle can be computed from the 
following formulation of convective heat transfer model (Hsu and 
Rubinsky, 1987) : 

Nu = 0.2233 (X + ) -0 1 745 5 + 3.66 (4-3) 


where 


X 1 

X + = — • 

R (Re • Pr) m 

Here 

X = Axial Distance, half size of nozzle length is used for 
calculation 
R = Radius of Nozzle 
Re = Reynolds Number of Plasma 
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Pr = Prandtl Number of Plasma 


Nu = Nusselt Number of Plasma 
Subscript m = mean values 

Heat convection loss within the nozzle, thus, can be computed from the 
Nusselt number obtained from Equation (4-3) , namely, 

^ (hp lasni “ hjl ) 

Q*,c " Nu • A„ • - • (4-4) 

D C P 

where 

Q m t c = Heat Convection Loss within the Nozzle 
A„ = Heat Transfer Area of Nozzle 
D = Diameter of Nozzle 

k = Thermal Conductivity of Plasma Arc 
Cp = Constant Pressure Specific Heat of Plasma Arc 
^plasma = Average Enthalpy of Plasma Arc through the Nozzle 
h M = Average Enthalpy of Nozzle corresponding to the Nozzle 

Wall Temperature 

All thermal properties of physical parameters used in this model shall 
be checked from the temperature corresponding to the average enthalpy 
of plasma arc. 

(IV-C-2) Direct Power Loss from the Surface of Nozzle 

Direct power loss from the nozzle surface includes P M / i , nozzle 
work function induced power input; and part of P H / 2 , nozzle anode and 
cathode drops induced power input, namely 

<2x,G = P*,l + ^X,2 (4-5) 

where , 

Q K/G = Direct power loss from the nozzle surface 
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i = Fraction of nozzle power loss contributed by anode and 
cathode drops 

(l-/)P K/2 = Plasma arc power input contributed by nozzle anode 
and cathode drops 

(XV— C— 3 ) Electron Collision Power Loss Within Nozzle 

Power loss conducted into nozzle by free electrons can be 
calculated from the following equation: 


Qx 


e 



e 


3 \ 

- KT 

12 J 


(4-6) 


where 

Qji e = Free electron power loss in the nozzle 
e = Electron charge = 1.602 x 10" 19 c 0 u 1 om b / m 0 1 e 

K = Boltzman constant = 1.38 x 10“ 16 erg / roo i*K 

T = Plasma temperature of electrons entering nozzle, 
approximately; and I p is the mean pilot current which can be computed 
as follows: 


t + + t_ 

Total power loss in nozzle, Q H , is the summation of Equations (4-4) , 
(4-5) , and (4-6) , namely, 

Qm = Qx,c + Qx,g + Qx,e ( 4 “ 8 ) 

( IV— D ) Power Loss in the Standoff Column 

Power loss in the standoff column includes radiative heat loss, 
Q S/r , and convective heat loss, Q S/C . 

Calculation of radiative heat loss is similar to Equation (4-2) 
for that in the gap between electrode and nozzle, namely 
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(4-9) 


Q S/ r( w )= V s (m 3 ) (4 x 10 10 


T(°K) 


11 6 


15000 


1 + 


T(°K) 


15000 


1 6 ' 


where V s denotes plasma arc volume within the standoff column; and T is 
the temperature corresponding to the average value of plasma arc 
enthalpy. 

( IV— E ) Total Workpiece Power Loss 

The shape of plasma contact areas in the keyhole of the workpiece 
is taken to resemble a finger ring. Computations yield a heat 
transfer rate out of rear half of a keyhole rather than that of the 
front half. This ring-shaped geometry can be simplified and expanded 
as illustrated in Figure 2. 

The keyhole surface has been divided into 5 equal segments along 
the keyhole axis between the inlet and outlet ends. Figure 2 shows 
the geometrical parameters of the segment. Effective heat transfer 
areas of the keyhole can be computed as follows: 

L = L 4 /5 
a t = irr * 


a = tan 


- i 


fl 

-L, 


Tt 

— r 

12 ) 


= tan 


- i 


,*r / 


L 

a, = a, + 2 

tan a 0 

1 

Aj = A s = — (a 1 + a 2 ) L 
2 


itr + 


L ' 
tan a°j 


L 


A 2 — A 4 — 


a t + 2 


L ) 


tan a 
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2L 


7vr + 


A 3 = 


tan a 0 ) 
L ' 


2tr - 


tan a 0 / 


where 

r = Radius of keyhole 

L 4 = Thickness of workpiece 

and A : , A 2 , A 3 , A 4 and A s denote effective heat transfer areas of 
keyhole at Segments 1, 2, 3, 4, and 5, respectively. 

Total power loss in the workpiece is comprised in the following 
items : 


(IV-E-1) Convective Heat Loss in Workpiece 

Governing equation is same as Equation (4-3) , used in nozzle, 
which is 

Nu = 0.2233 (X + ) "° • 74 5 5 + 3.66 (4-10) 

where 


x 1 

x + = 

R (Re • Pr) m 

Here x is the half thickness of workpiece to be used as the distance for 
calculation, and the rest of symbols stand for the same meaning as that 
used in the nozzle. 

Heat convection loss in the workpiece can be computed from the 
Nusselt number obtained from Equation (4-10) , namely, 

^ (^plasma ~ hy) 

Q H/C = NU • A v • (4-11) 

D C P 

where A y stands for heat transfer area of workpiece which is the shape 
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of plasma contact area in the keyhole and is taken to resemble a finger 
ring (see Figure 2 for details) ; and D is the diameter of keyholes 
which is also the diameter of nozzle if plasma arc jet disperse through 
standoff column is small enough to be neglected. The rest of symbols 
are the same as that used in the nozzle. 

(IV-E-2) Direct Power Loss from the Surface of Keyhole 

Direct power loss from the keyhole surface includes ;? W/ . , 
workpiece work function induced power input; and part of ? V/J . 
workpiece anode and cathode drops induced power input, namely 

Cv,G = Pv,! + >»Pv , 2 < 4 - 12 ) 

where 

Q W/G = Direct power loss from the nozzle surface 

+ = Fraction of workpiece power loss contributed by anode and 

cathode drops 

( 1— ^ ) Py 2 ~ Plasma arc power input contributed by workpiece anode 
and cathode drops 

In the present model, it is assumed that ^ = 0. This is based on 
the assumption that both anode and cathode drops occur in the plasma 
flow. Power loss generated by anode and cathode drops are thus 
absorbed by the plasma flow which raises up the average enthalpy of 
plasma jet. This higher value of plasma jet enthalpy is the 4 -, 
transferred to the workpiece through the convective heat transfer. 

In this calculation, the plasma flow parameters used in the model 
are all one-dimensional. This means that the plasma flow parameters , 
such as power loss, enthalpy and temperature vary with respect to th:.i 
axial coordinate of flow direction, while the average values are taken 
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for all parameters transverse to the flow direction. This assumption 
matches the concepts described earlier. 

(IV— E— 3) Electron Collision Power Loss in Keyhole 

Power loss due to free electron collision on the surface of 
keyhole can be calculated from the following equation: 

f 3 ) 


- KT 


(4-13) 


e [2 ) 

where Q V/ e denotes free electron power loss in keyhole; T, the average 
temperature between plasma arc and wall of keyhole; and I m , the 
electric current which can be computed from the following equation: 
Im + t + 


(4-14) 


t+ + t. 

Thus, total power loss in workpiece, Q w , is the summation of Equations 
(4-11), (4-12), and (4-13), namely 

Qw - Qv,c + Qy,0 + Qv,e (4-15) 

(IV-F) Power Discharge Loss 

Amount of power loss carried away by the plasma jet discharge 
from the exit end of keyhole, Q d j s , can be computed from the plasma arc 
enthalpy, h exit , multiplied by the mass flow rate of argon gas, m, 
namely 


Qd i s x i t * m (4 — 16) 

Obviously, summation of power loss throughout the entire VPPA 
welding system is equivalent to the total power input of the system, 
p t o t a l ' namely 


P total = Qe + Qg + Qx + Qs + Qw + Qd i s (4-17) 

Calculation of power loss at each part of VPPA welding system 
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will be illustrated in Section VIII of this report. 
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V. Calculation of Plasma Arc Enthalpy 


Density of argon gas at the pressure of 1 atmosphere and 

temperature of 298 °K is 

PM Ar 101320 (N/m 2 ) • 39.95 (kg/kmol) 

P ~ = 

Ar R 0 T 8315 (N-m/kmol - °K) • 298 (°K) 

= 1.6336 (kg/m 3 ) 

where 

P = Atmospheric pressure 
M A r = Molecular weight of argon 
R 0 = Universal gas constant 
T = Temperature 

Mass flow rate of argon gas in terms of volumetric flow rate can be 
shown as 


m = (V ft 3 /hr) 


r 1 hi ' 

hr/s m 3 /ft 3 


13600 


135.3147 


A r 


1.6336 


(3600) (35.3147) 


(V ft 3 /hr) 


1.285 • 10" 5 (V ft 3 /hr) (kg/s) 


1.285 • 10" 2 (V ft 3 /hr) (g/s) 


A shielding gas is needed to isolate the welding system from the 
surrounding environment. Shielding gas purity is an essential factor 
in achieving quality welds in any gas shielded arc welding process. 
The problem of impurities, when manifest, is normally resolved by 
eliminating the source of the contamination. However, if the gas 
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supply itself is contaminated, serious schedule interruptions can 
result, particularly where large manifold gas distribution systems 
are employed. 

Gas contamination affects VPPA welding in several ways. Inside 
the torch (standoff column) , oxygen from water vapor decomposed in the 
arc can produce a visible oxide film stain on the tungsten electrode. 
The plasma jet may sputter and exhibit a green coloration from copper 
eroded from the nozzle under erratic flow conditions. The weld 
itself may exhibit a packed-looking surface, presumed to be fine 
porosity due to hydrogen gas from decomposed water vapor. Under 
severe conditions, the torch may leave behind it a train of 
irregularities and holes in the weld bead. 

Plasma arc power distribution in the whole VPPA welding system 
can be divided into the following five parts for convenience of 
discussion: (a) Plasma arc power within the gap between electrode 
and nozzle, (b) plasma arc power within nozzle, (c) plasma arc power 
within standoff column, (c) plasma arc power within keyhole of 
workpiece, and (d) plasma arc power discharged from the workpiece. 
The plasma arc power associated with ionized argon gases is 
characterized by the thermal energy parameter enthalpy, which is a 
function of temperature. The relationship between enthalpy ancl 
temperature is non-linear for high temperatures as encountered in the 
plasma. For these computations the plasma enthalpy-temperature 
relation is obtained from the tables for argon computed by Drellishak 
et al. (1963). Thermal properties and transport coefficients of 
argon plasma are also a function of temperature. The high 
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temperature argon gas transport coefficients can be obtained from 
Lancaster (1986), Metcalfe and Quigley (1975), etc. As to aluminum 
as the workpiece for the present study, the thermal properties can be 
obtained from American Society for Metals Handbook Committee (1979) . 

The itemized calculation of plasma arc enthalpy at various 
sections of VPPA welding system are illustrated as follows: 

(V-A) Plasma Arc Enthalpy Within the Gap Between Electrode and 
Nozzle 

Let us assume that entrance and exit ends of plasma arc enthalpy 
in the gap between electrode and nozzle are h 1;i and h 1/0 , 
respectively . 

The plasma arc power at the entrance end of the gap, P x i ', is 
power input at electrode, P E , subtracting electrode heat conduction 
loss, Q e , namely 

p i,i = p e - Qe (5-1) 

Plasma arc enthalpy at entrance end of the gap, thus, can be calculated 
as 



m 


The plasma arc power at the exit end of the gap, P 1/0 , absorbs 
power through Joule heating, Pj , l , and gives away power through 
radiative heat loss, Q G . It gives 

p i , o = p i,i + p j • i “ Qg (5-3) 

Thus, plasma arc enthalpy at exit end of the gap can be calculated from 
the expression 
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(5-4) 


hi 


0 


Pi 


0 


m 


(V-B) Plasma Arc Enthalpy Within Nozzle 

Let us assume that the plasma enthalpy at the entrance and exit of 
the nozzle are h 2/i and h 2/0 , respectively. 

The plasma arc power, P 2 / j , and enthalpy, h 2/i , at the entrance 
the nozzle are identical to the plasma arc power, P 1/0 and enthalpy, 

hi , o the exit end of the gap between electrode and nozzla, 
respectively. Thus 


Plasma arc power at the exit end of nozzle, P 2 / 0 , absor. . power 
from partial contribution made by nozzle anode drop, (l-^)P N 2 , and 
full contribution made by Joule heating within nozzle, P, , , and 
subtracts convective heat loss in nozzle, Q„ # e , and free electron 
power loss in nozzle, Q„ . Thus, 

P 2,o = P 2 ,i + (W)P» /2 + Pj, 2 * Qx,c - Qx, e (5-6) 
The corresponding plasma arc enthalpy at the exit end of nozzle is 


2 , 0 


2 , o 


m 


(5-7) 


(V-C) Plasma Arc Enthalpy Within Standoff Column 


Let us assume that entrance and exit ends of plasma arc enthalpy 
in the standoff column are h 3 ^ and h^ / 0 / respectively. 

The plasma arc power, P 3 / j , and enthalpy, h 3 f j , at the entrance 
end of standoff column are identical to the plasma arc power, P ? 0 , rind, 
enthalpy, h 2/0 , at the exit end of nozzle, respectively. Thus, 
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(5-8) 

^3 / i — ^2 , o ' 

Plasma arc power at the exit end of standoff column, P 3 7 0 , 
absorbs Joule heating in standoff column, P J/3 , and subtracts 
radiative heat transfer loss in the same section, Q s . Thus, 

P 3 , o — P 3 , i + P J , 3 ” (5-9) 

The corresponding plasma arc enthalpy at the exit end of standoff 
column is 


h 


4 , 0 



(5-10) 


m 

(V-D) Plasma Arc Enthalpy Within Keyhole of Workpiece 

Let us assume that entrance and exit ends of plasma arc enthalpy 
in the keyhole of workpiece are h 4 t ^ and h 4 t 0 , respectively. 

The plasma arc power, P 4 t j , and enthalpy, h 4 t j , at the entrance 
end of the keyhole of workpiece are identical to the plasma arc power, 
P 3 1 0 , and enthalpy, h 3 / 0 , at the exit end of standoff column, 
respectively. Thus, 

P 4 , i = p 3 , o 

^4 , i — ^3 , o 

Plasma arc power at the exit end of the keyhole of workpiece, 
P 4 ' Q , absorbs power from partial contribution made by workpiece anode 
and cathode drops, (l-/)P U/2 , and subtracts convective heat loss in 
keyhole, Q W/C and free electron power loss in keyhole, Q W/6 . Thus, 

P 4,o = P 4,i + ( 1”/ ) p w , 2 - Qv,c - Qv,e ( 5-12 ) 

The corresponding plasma arc enthalpy at the exit end of the 
keyhole of workpiece is 


(5-11) 
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( 5 - 13 ) 


(V-E) Plasma Arc Enthalpy Discharged from the Workpiece 

Plasma arc power, Q dis , and enthalpy, h d { s discharged from the 
keyhole of workpiece is identical to the plasma arc power, P 4/0 , and 
enthalpy, h 4/0 , at the exit end of keyhole, namely 


Qd is ^4 , 
^d i S = ^4 , 


( 5 - 14 ) 


Calculation of plasma arc enthalpy at each part of VPPA welding 
system will be illustrated in Section VIII of this report. 
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VI. Calculation of Crown and Root Bead Widths 


Bead width is calculated from the heat absorbed in the workpiece, 
Q h . As shown in Equation (4-15) Q M is the summation of convection 
heating (The empirical relation actually includes radiation.) from 
the plasma, "direct" heating of the surface by current passing through 
localized surface voltage drops, and heat carried into the surface by 
hot electrons from the plasma. 

For the purpose of estimating the melting isotherm surface of the 
weld puddle the total heat absorption Q w is taken to enter the 
workpiece at a distance $L 4 above the bottom of the plate, where L ; 
stands for the thickness of the workpiece, and (j) represents a fraction 
of the thickness. If $ = 0.5, the heat source is located at the cents:-, 
of the plate and the temperature distribution is symmetrical such that 
the puddle crown and root will be identical in size: an appropriate 

model for a full penetration keyholeing plasma arc weld. If <J) >0.5, 
the heat source is located closer to the crown side of the plate and the 
crown will be wider than the root: a better model of a typical 

keyholeing weld. If <J> = 1.0, all the heat is delivered at the crown 
surface: a model suitable for representation of a GTA weld or a 

partial penetration VPPA weld. In the illustrative computations 
included in this report (J) = 0.55. By adjustment of ({) it is possible to 
model effects of real assymetrical distributions of heating in the 
keyhole column as well as the further redistributions of heat due to 
circulations in the weld puddle and thus to model approximately the 
entire range of VPPA penetration phenomena. 
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The temperature field within the confined space between the two 
surfaces (assumed nonconducting) of the workpiece is constructed out 
of moving heat sources in an infinite continuum. ''Image" sources are 
added outside the workpiece space so as to produce planes of symmetry 
across which no heat flows at the location of the workpiece surfaces. 
Given coordinate z zeroed at the bottom of the workpiece and extending 
up into the workpiece body and locating the heat source Q„ at z = <|)L 4 , 
then an image heat source below the bottom of the plate of strength Q H 
located at z = <J)L 4 creates a plane of zero heat flow at z = 0, the 
location of the bottom surface of the workpiece. The latter plane 
will remain and a further plane of zero heat flow at the site of the top 
surface of the workpiece will be created if an infinite array of image 
heat source pairs is constructed along the z-axis, each heat source 
pair of strength Q H spaced at ±<J)L 4 about an infinite array of points 
spaced at 2L 4 and including z = 0. 

The temperature T c at location (x c , y c , z c ) within the workpiece 
due to the welding power input to the workpiece is then: 

. Qh f u(R| + x c n 

T c - T 0 =. E exp (6-1) 

1 = 0 4 7TKR j l 2a 

where 

T 0 = Ambient temperature of workpiece 
K = Thermal conductivity of workpiece 
u = Weld speed (in the x-direction) 
a = Thermal diffusivity of workpiece 

R, = Distance between the i th point heat source and the point 
of interest 
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Rj is computed as follows: 


Ri = -J x c 2 + y c 2 + (z t - z c ) 2 ( 6 - 2 ) 

Weld bead width is determined by locating the largest value of y c 
for which T c = T^,, the melting temperature of the workpiece, at z = 0 
(root) and z c = L 4 (crown) . The weld widths are twice the values of 
these y e 's. A converging computer search routine is used to make the 
computation. 

Two approximations are used in the computations. First, the 
summation is truncated. 3 5 source points are used, 17 image sources 
on each side of the weld heat source. The values of z - x are generated 
by a recurrence relation: 


Z j — Z i _ 2 2 L 4 

(6-3) 

where the starting values are: 


2 ^ = -I 6 L 4 (J) L 4 

(6-4) 

~~ t (j)L«4 

(6-5) 


An evaluation of the Zj's will show that z 18 (=<j)L 4 ) is the weld heat 
source located in the workpiece. 

Second, although the maximum values of y c on the melting isotherm 
occur somewhat behind x c = 0 , for the sake of a substantially easier 
compuation it is considered sufficiently accurate to accept for the 
maximum value of y c the slightly (at the values of u appropriate to 
VPPA welding) smaller value of y c that intersects the melting isotherm 
at x c = 0 (Solomon and Levy, 1981; Tsai, 1982) satisfying the equation: 



1 = 0 4 ttKR j 


exp 



(6-6) 


where 
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R i ■ "J Yo 2 + ( 2 i " z c) 2 

z e = 0 for the root 

z e = L 4 for the crown 

and the are generated in accordance with Equations (6 
(6-5) . 


(6-7) 
(6-8) 
(6-9) 
3) to 
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VII. Calculation of Crown and Root Heights 


(VII— A) Volume Conservation Equation of Materials 

Reinforcement profiles of crown and root beads are represented 
as paraboloids, as shown in Figure 4. Paraboloidal profiles can be 
expressed as a special case (where n=2) of the following equation: 

(7-1) 




'2y' 

n “ 

Z = h 

1 - 





id j 



where 


y = Distance from weld centerline 
z = Local bead height 
h = Maximum bead height 
d = Bead width 

n = Numerical value greater than 1 
The cross-sectional area of the crown reinforcement is 


A c 



r 2y> 

n ■ 

IdJ 



dy 


n 

= h c d c 

n + 1 

where 


h c 

<*c 

Similarly, 

A r 


= Crown height 
= Crown bead width 

cross-sectional area of the root reinforcement is 
n 

h rd r 

n + 1 
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where 


h r = Root height 
d r = Root bead width 

Volume conservation equation can be established based on the 
equilibrium relation between the material volume provided by welding 
wire and reinforcement volume enlargement of crown and root beads 
which is given by 

n n n V y 

d c h c + d r h r « - d v 2 L 4 • W g (7-2) 

n+l n+1 4 u 


where 


d w = Diameter of wire 
V w = Wire feed rate 
u = Weld speed 
W g = Initial gap at weld seam 
L 4 = Thickness of workpiece 
(VII— B) Force Equilibrium Equation 

For the purpose of deriving the force equilibrium equation, one 
has to calculate cross-sectional areas a, b, and c, defined in Figure 
5, which are given by 


A a 


Similarly, 

A b 




h c d r - 


n+l 


h c d c 


fd r ) n+1 


\U,/ 


d r • L 4 


(7-3) 


(7-4) 
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and 


n 


A c = A r = 


h r d r 


(7-5) 


n + 1 

Slopes at points B (crown) and D (root), shown in Figure 5, can 
be expressed as follows: 

2nh c 

y=d r /2 ' d c 


dz 

tan e c = — 
dy 


d r ) n_ 1 


and 


tan 0 r = 


dz 


dy 


2nh, 


y=d r /2 d r 

With reference to the free-body diagram shown in Figure 5 , one jan 
derive the force equilibrium equation by the balancing of forces, 
which include gravitation, surface tension and pressure acting on the 
molten weld puddle: 


' 

1 


[ d r ] 

i — 

G 

H 

+ 

c 

h c d r - 

h c d c 



+ d r • L 4 + d r h r 

■ 

n + 1 


id c i 


n + 1 


-h r h c 

[<V 

n - : 

L " 


= 4n 

— 




<r - p d r 


Ur d c 






p g 


(7-6) 


(VII -C) Calculation of Crown and Root Heights in Flat 
Welding (with Consideration of Gravity Force) 

In comparison between the effects of external pressure p and 
surface tension force acting on the molten fluids of welding, 
experience suggests that the external pressure is much smaller than 
that of the surface tension force. Therefore, p will be ignored in 
the force equilibrium equation (7-6) . Crown and root heights, thus, 
can be solved simultaneously from Equations (7-2) and (7-6) which are 
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given by the following expressions: 


4 <r 

1 2 (n + 1)~ 


h c = 


dr' 


P<3 


It V„ 

- d v 2 L 4 

4 U 


W, 


- d r -L 4 


n 


d r " 


n + 1 


d c + 


4n<r 

' i 

f d r] 

d c ' 
^ 1 " 

1 

f d r| 

p g 

. d c 

^ d c j 

-3r J J 

a c 

n + 1 

id c j 


n + 1 


(7-7) 


tv V v n 

- d v 2 L 4 ■ Wg 

4 u n + 1 


d c h c 


h r = 


(7-8) 


n 


n + 1 

(VII-D) Calculation of Crown and Root Heights in Vertical Welding 
(Without Gravity Effect) 

The effect of gravity force on the molten fluids of welding can be 
ignored for the case of vertical welding or the welding carried out at 
in space at zero-gravity. In this respect, one can ignore the terms 
of both gravity force and external pressure in Equation (6-6) which is 
given by 


h r = h c 


(d T \ n 


'd c 


(7-9) 


Substituting this equation to Equation (7-2) , on obtains 
TV V„ 


- <V L, • W, 

4 U 


h c = 


(7-10) 


n 


n + 1 


d c + d r 


ld c J 


(VI -E) Solidification of Crown Height 

Figure 6(A) shows the configuration of the cross-section of a 
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puddle with the interface between liquid (molten metal) and solid 
metal. As solidification occurs interfaces between the molten and 
solid metals of the workpiece move together. Figure 6(B) shows the 
configuration of puddle with a trapezoid liquid-solid interfaces in 
contact at the just-completely solidified puddle root. It is 
postulated that the root solidifies before the crown because the root 
is farther from the heat source and, usually, the root is narrower. 
Figure 6 (C) shows the configuration of the cross-section of the puddle 
after complete solidification. This figure shows shrinkage of the 
crown height that took place during solidification of the reservoir of 
molten metal above the frozen root. The cross-sectional area of the 
residual molten metal just after root solidification can be 
approximated as 
1 

A c = — L 4 • d c ' (7-11) 

2 


where L 4 stands for the thickness of the workpiece; and d c ' , the crown 
width of the remaining reservoir of molten metal. Given specific 
volumetric shrinkage, av/v, of the metal during solidification a 
change of cross-sectional area of the remaining puddle molten zone 
takes place according to the followin relation: 


AA c = A, 


' AV' 

1 

' AV 1 


- - L 4 • d e ' 


lv J 

2 



(7-12) 


From Equation (7-5) for the cross-sectional area of crown 
reinforcement, variation area of crown reinforcement can be shown as 


A A c = A 


n 


\n + 1 


<Vhc' 


n 


n + 1 


d c 'Ah c ' 


(7-13) 


41 


Equating Equations (7-12) and (7-13) leads to 


- Li 


<V 


n 


V. V t 


n + 1 


d c 'ah c ' 


or, simplifying, to 


n + 1 ( av) 

Ah c ' = L 


2n 


\v 


(7-14) 


According to King (1987) , the specific volumetric shrinkage for 
aluminum during the solidification, av/v, is 6.7%. Therefore, crown 
height shrinkage for aluminum during solidification can be expressed 


as 


n + 1 

Ah c ' = 0.067 L 4 (7-15) 

2n 

Detailed calculation results of crown and root heights are 
illustrated in Section VIII of this report. 
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VIII. Results of Computation 


Recently, some limited experiments relating VPPA welding system 
parameters to voltage and bead configuration outputs have been 
carried out by General Digital Industries (GDI) through the 
sponsorship of NASA's Marshall Space Flight Center. Some of the 
input parameters used by GDI have been adopted in this study to permit 
comparison of the experimental observations made by GDI and our model 
computations. Tables 1 and 2 contain the basic input parameters to be 
adopted in the model computation. By using these input parameters, 
all the physical parameters listed for the VPPA welding system are 
computed . 

Table 3 shows the electric potential illustrated in Figure 1 for 
standoff distances from 1 mm to 10 mm as computed by the model. All 
the symbols shown in the table are illustrated in the attached table of 
nomenclature . 

Table 4 shows the itemized power inputs* at various sections of 
VPPA welding system for standoff distances from 1 mm to 10 mm from 
model computation. 

Table 5 shows plasma jet power, enthalpy and temperature 
profiles at various locations of VPPA welding system for standoff 
distances from 1 mm to 10 mm from model computations. 

Table 6 shows power loss at various sections of the VPPA welding 
system for standoff distances from 1 mm to 10 mm from model 
computations. From this table we draw the following conclusions: 

(a) Total electric power consumption increases as the standoff 
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distance increases; 

(b) Power discharged from the outlet end of keyhole increases as 
the standoff distance increases while the percentage of the 
ratio of power discharged from the outlet end of the keyhole 
to total power consumption decreases as the standoff 
distance increases; 

(c) Power absorbed by the workpiece, which melts metal and 
produces the weld, increases as the standoff distance 
increases as does the power discharged from the outlet end of 
keyhole ; 

(d) Power absorbed by the workpiece is less than the power 
discharged through the outlet end of keyhole; 

(e) Power loss percentage in the gap between electrode and nozzle 
and at the nozzle both decrease as the standoff distance 
increases while the power loss percentage in the standoff 
column increases ; and 

(f) Plasma jet temperature at the outlet end of the keyhole 
ranges from 13700 to 14900 °K as the standoff distance varies 
from 1 mm to 10 mm. 

Table 7 shows crown width, root width and average width of 
workpiece computed from the model computation. 

In the GDI experimental data available to date welding wire 
effects were not yet included. Hence these effects were not included 
in our calculations. 

Computed results of total power input, power discharged, 
workpiece power loss and standoff power loss for various standoff 
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distances from 1 mm to 10 mm are shown in Figure 7. Figure 8 shows 
computed results of percentages of standoff jet column power loss, 
workpiece power loss and plasma discharge power loss for standoff 
distances from 1 mm to 10 mm. 

Comparisons between the results of model computations and 
experimental observations are shown from Figures 9 to 13. Figure 9 
shows computed transferred arc voltage at straight polarity r.c: 
standoff distances from 1 mm to 10 mm with experimental observations 
for comparison. 

Figure 10 shows computed transferred arc voltage at reverse 
polarity for standoff distances from 1 to 10 mm, with experimental 
observations for comparison. Figures 11, 12, and 13 show computed 
crown width, root width, and average bead width, respectively, fo 
standoff distances from 1 to 10 mm, in comparison with that, of the 
experimental observation. These results shown in Figures 9 to 13 
clearly indicate that the computed results are in excellent agreement 
with the experimental measurements. 
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IX. DISCUSSION AND CONCLUSION 


The purpose of this study is to provide an algorithm for 
computing crown and root widths and crown and root heights of a VPPA 
weld bead to be used as a component in a VPPA welding control system. 

The power distribution has been analyzed for the argon plasma gas 
working fluid in the VPPA welding process. 

The input parameters used to check this model are listed in 
Tables 1 and 2 . Some of the input parameters were provided by GDI . 

In the GDI experimental data available to date welding wire 
effects were not yet included. Hence these effects were not included 
in the calculations. Figures 9 to 13 show clearly that computed 
results of transferred arc voltages at both straight and reverse 
polarities, and crown, root and average bead widths for standoff 
distances from 1 to 10 mm are so far in good agreement with the 
experimental measurements. 

It must be remembered that the model as it stands now is not 
complete. For example, it remains to incorporate into the model the 
effect of surface tension gradient driven circulations, which are 
extremely sensitive to surface contamination levels of various sorts. 
Small differences in chemistry can result in great differences in 
puddle configuration and weld penetration (Lambert, 1991) . 

The most important and critical points for the model 
computations are that some physical parameters and their profiles, 
which are not measured by the experimental observations, can be shown 
and illustrated clearly from the mathematical computations. These 
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important physical parameters and their profiles include power input, 
power loss, enthalpy and temperature profiles and distributions at 
various sections of VPPA welding system. Also, what are the main 
causes and contributors to provide the mechanisms of power loss can 
only be disclosed from the model computation. These important 
physical mechanisms are illustrated in Tables 4 to 6 and Figures 7 and 
8 . 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2 . 

Figure 3 . 

Figure 4 . 
Figure 5 . 


Figure 6. 


Figure 7 . 


Figure 8 . 


Figure 9 . 


VPPA electric power system. The left side shows the 
straight polarity mode, and the right shows the reverse 
polarity mode. 

Geometry and dimensions of the keyhole's approximate, 
simplified heat transfer area. 

Point heat source in workpiece and 17 image heat sources 
added above and below. 

Dimensions of the weld puddle cross-section. 

Free body diagram for the balance among external 
pressure P, surface tension force <r, and gravitational 
force g over an approximate weld puddle cross section. 
Conjectured configuration of the puddle cross-section 
during solidification. (A) Fully molten cross-section 
before phase transition, (B) Partially molten cross- 
section just after freezing of the root; and (C) 
Completely solidified cross-section. 

Plasma arc total power input, power discharged, 
workpiece power loss and standoff power loss in terms of 
various standoff distances from 1 to 10 mm. 
Percentage of standoff power loss, workpiece power loss 
and plasma arc jet discharged power loss in terms of 
various standoff distances from 1 to 10 mm. 

Computed and observed main arc straight polarity 
voltages in terms of various standoff distances from 1 
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to 10 mm. 


Figure 10. 


Figure 11. 


Figure 12 . 


Figure 13 . 


Computed and observed main arc reverse polarity voltages 
in terms of various standoff distances from 1 to 10 mm. 
Computed and observed crown widths in terms of various 
standoff distances from 1 to 10 mm. 

Computed and observed root widths in terms of various 
standoff distances from 1 to 10 mm. 

Computed and observed average bead widths in terms of 
various standoff distances from 1 to 10 mm. 
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Table 1 Basic Input Parameters (I) 


Description 


Data Input 



Distance between electrode and nozzle 


Vertical length of nozzle 


Thickness of workpiece 


1.143 mm 


3.378 mm 


9.525 mm 


Length of electrode 


Nozzle diameter 


3.175 mm 


Electrode diameter 


3.9688 mm 



Thermal conductivity of electrode 


Argon gas mass flow rate 


Electrode work function 


0.1028 g/s 


2.6 V 


Nozzle work function 


4.5 V 


Workpiece work function 


4.2 V 


Anode potential drop 


4.4 V 



Cathode potential drop 


Main electric current in straight polarity 


8.8 V 


200 A 


Main electric current in reverse polarity 200 A 


Pilot electric current in straight polarity 


Pilot electric current in reverse polarity 


Time duration for straight polarity 19 ms 


Time duration for reverse polarity 



Voltage drop per unit length 


2 V/mm 



Electrode temperature 


Nozzle temperature 


Environment temperature; Initial 
temperature of workpiece 


2500 °K 


916 °K 


298 °K 

























































Table 2 Basic Input Panneters (II) 


Description 


Data Input 





Plasma enthalpy corresponding to nozzle 
wall temperature 


Plasma enthalpy corresponding to keyhole 
temperature 


Fraction of power loss due to anode and 
cathode drops 


Location of welding heat source from 
back of workpiece in fraction of plate 
thickness 






Welding travel speed 


Surface tension 


Feed wire diameter 


Wire feed speed 






Initial gap distance between weld surfaces 
to be joined 


Thermal conductivity of workpiece 


Specific heat of workpiece 


Density of workpiece 


Exponent of reinforcement shape 


kJ 

0.3672 — 

g 


kJ 

0.3672 

g 



N 

0.865 - 


1.5748 mm 


0.000 


0.000 mm 


0.13 


0.864 


m- 0 K 


2650 kg/m- 

































Table 3 


Computed Electric Potential 


Standoff (m 

u) i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

v E 

(V) 

-8.6 


-8 . 6 

-8.6 

-8.6 

-8.6 

-8.6 

-8 . 6 

-8.6 

-8 . 6 j| 

v c 

(V) 

-17. 3« 

SB 

IMP 

ifettl 

B 

■0 

■Si 

B 

B 

mm 

g 



(V) 

-10.6 


-14 . 6 

g 

-18.6 

-20.6 

-22.6 


g 

m 

V B 

(V) 



BB 

B 

■m 

mm 

-29.64 

g 

iSb 

B 

1: 

1-37.64 j 

)• 

j< 

V D 

(V) 

-8 . 46 

Si 

-12 .46 

a 

-16.46 

-18.46 



-24 .46 

5—2 6 . 46 f 

v A 

(V) 

-25.84 

-27 . 84 

-29.84 

-31.84 

-33.84 

-35.84 

-37.84 

-39.84 

-41.34 

-43 . 84 

' 

■j AV n. + 

(V) 

25.84 

27.84 

29 . 84 

31.84 

33.84 

35.84 

37.84 

39.84 

41.84 

43.84 | 

l! 

aV p + 

(V) 

17.39 

17.39 

17.39 

17 . 39 

17.39 

17.39 

17.39 

17.39 

17.39 

1 

17.39 

v E ' 

(V) ! 

16.08 

16.08 

16.08 , 

16.08 

16.08 

16.08 

16.08 

16.08 

16.08 

16.08 j- 

| 

I 

*c' 

(V) : 

24.84 

26.84 

28 . 84 

30.84 

32.84 

34.84 

36.84 

38.84 

40.84 

42.84 | 

1 

V 0 ' 

(V) 

18 . 08 

20.08 

22 . 08 

24.08 

26 . 08 

1 

f 

28.08 

30.08 

32 . 08 

34 . 08 

36.08 ! 

i 

j 


(V) ; 

n 



B 



39.12 

41.12 

43 . 12 

45.12 

v d' 

(V) ; 

26.98 

28.98 

30.98 

32.98 

34.98 


38.98 

40 . 98 

42.98 

44.98 

v a' 

(V) 

34.12 | 

36.12 

38.12 | 

40.12 

42 . 12 

44.12 

46.12 j 

48 . 12 

50.12 

52 . 12 

pv m * 

(V) 

34 . 12 

36 . 12 

38.12 j 

40 . 12 

42 . 12 

44.12 - 

46.12 

48.12 

50 . 12 

52.12 

1 

Ou 

> 

<1 

(v) 1 

1 

-7.14 

-7.14 

-7 . 14 .■ 

-7.14 j 

i 

-7.14 j- 

-7.14 • 

-7.14 j' 

-7.14 • 

-7 . 14 
— 

-7 . 14 










































Table 4 


Computed Power Input 



Standoff (nun) 


9 10 



1341.1 


P N,1 

(W) 

P *,2 

(W) 

P W,1 

(W) 

P W , 2 

(W) 


128.2 

128.2 

128.2 

128.2 


125.3 125.3 125.3 125.3 125.3 


(W) 1033 1033 jl03 3 1033 1033 1033 -■ 1033 1033 1033 1033 


486.5 1*86.5 i486. 5 486.5 486.5 486.5 486.5 486.5 436.5 
! 1 - 


1298 . 41298. 4|l298. 4 1298 .4 


Pj , t (W) 486.5 486 


1298 . 4 


400 800 [1200 h.600 2000 2400 2800 3200 3600 4000 


Total 


8559 . 68959 . 


II 


































Table 5 Computed Plasma Arc Jet Power, Enthalpy and Temperature 


Standoff (mm) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

P l,i (W) 

1244 . 1 

1244.1 

1244.1 

1244.1 

1244.1 

1244 . 1 

1244.1 

.1244.3 

.1244 . 3 

.1244.1 

h t/i (kJ/g) 

12 . 1 

12.1 

12.1 

12.1 

12.1 

12.1 

12 . 1 

12 . 1 

12 . 1 

12 . 1 

*1,1 (°K) 

12046 

12046 

12046 

12046 

12046 

12046 

12046 

12046 

12046 

12046 

p i,o'-P 2/ i (W) 

1713.5 

1713 . 5 

1713 . 5 

1 

1713 .51713 .5 

1713 . 5 

1713.5 

1713 . 5 

1713 . 5 

1713 .5 

k X , a ' k 2 , i ( kJ/ g ) 

16.67 

16.67 *16 .67 

16.67 jl6-67 

16.67 |l6.67 

16. 67 

16.67 

16 .67 

*l,o'**2,i (°K) 

12813 !].28 13 '12813 
1 1 

12813 

1 

12813 






' P 2,o' P 3,i (W) 

2389. 6|2389. 612389. 6 

2389. e! 

i 

2389.6 

2389.612389.6 

1 

2389 . 6 

2389 . 6 

2389. 6j 

mnn 

1 i 

23 .25 23.25 L 

i f 

23.25 j 

23.25 23.25 

1 

23.25 ^23.25 

i 

23 . 25 

2 3.25 

23 . 25 

T 2 /0 ;T 3/i (°k) 

i j 

13655 3.3 655 

i 

i I 

i ! 

13655 !13 655 13655 
1 1 





P 3 , o ' P 4 , i (W) 

2720.3!: 

i 

1 

JO27.4I: 

i 

3310.3 

3570.4: 

3807.6 

4024.54222.4 

B 

B 

4727.8 


26.46 29.45 *32-20 34.73 37.04 

; ; 1 

39.15 41.07 

42.84 


45.99 

BOD 

14035 14353 14646 14915 15182 

; ! : I 

1 1 > 1 

15435 15667 

! 

15880 

16115 | 

1 

i 

16375 j 

P 4,o (W) 

2458.8,2660. ll; 

i i 

• 1 

2837. 2)2 996. 63134. 2 

i 1 

J 

3254. 3j: 

3351.7 

■ 


h 

3559 .2j 

h 4/0 (kJ/g); 

1 j , 1 I ; 

23.92 25.83 ! 27.60 29.15 30.49 31.66 32.60 : 

'i'i!; 

33.41 

34.07 j 

34.62 [ 

t 

1 

— 4,o (°K) 

Ij ; 1 : 

13735 =13970 114156 14321 14464 14588 14689 1 

L4775 . 

14844 : 

L4903 } 
1 




























Table 7 Computed Bead Widths 



























Table 6 Computed Power Loss 


* 


2 

3 

97 . 0 

97.0 


Standoff (mm; 12 34567 89 10 


Q e (W) 97.0 97.0 97.0 97.0 97.0 97.0 97.0 97.0 97.0 97.0 


Qe/^T o t a 1 1-7% 



1.5% 1.4% 


bbbbb 


17.2 17.2 17.2 17.2 17.2 17.2 17.2 17.2 17.2 17.2 


Qg (W) 


Qg/Qto t a 1 0.3% 


Q„ (W) 


0 t a 1 


s (W) 69.3 162.2 279.1 



0.3% 0.3% 0.3% 0.2% 0.2% 0.2% 


875.7 


15 . 2% 



13 . 3% 12 . 6 


Qs/Qto t a 1 


Q H (W) 


Qh/Qt 0 t a 1 


i s 




2.6% 4.3% 6.0% 7.9% 9.9% 11 . 9% 13 . 9% 15 . 8% 17 . 8% 


2347. 4^2 453. 32553. 92 653. 


mm 


38.9% 38.1% 37.4% 






2996. 6(3 134. 23254. 3i3351. 713435 . 0i3 502 .13 559.2 


5.4% 

34.9% 

34.5% 

34.1: 






































































FIGURE CAPTIONS 


Figure 1. 


Figure 2 . 


Figure 3 . 


Figure 4 . 
Figure 5 . 


Figure 6 . 


Figure 7 . 


Figure 8 . 


VP PA electric power system. The left side shows the 
straight polarity mode, and the right shows reverse 
polarity mode. 

Ring-shaped geometry and dimensions of the keyhole's 
simplified heat transfer area. 

Distribution of real point heat sources in workpiece and 
17 points each of heat sources added to the mirror 
reflection corresponding image point sources at upper 
and lower planes 

Geometric coordinate of the weld puddle cross-section. 
Free body diagram and the geometry of coordinate to 
describe the balance among external pressure P, surface 
tension force a, and gravitational force g. 
Configuration of the puddle cross-section with the 
liquid (molten metal) -solid metal. (A) Molten state 
puddle before phase transition, (B) Molten state puddle 
in the process of phase transition from molten metal to 
solidified metal; and (C) solidified metal of welding 
puddle . 

Plasma arc total power input, power discharged, 
workpiece power loss and standoff power loss in terms of 
various standoff distances from 1 to 10 mm. 
Percentage of standoff power loss, workpiece power loss 
and plasma arc jet discharged power loss in terms of 
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Figure 9 . 


Figure 10. 


Figure 11. 


Figure 12 . 


Figure 13 . 


various standoff distances from 1 to 10 mm. 

Computed and observed main arc straight polarity 
electric voltages in terms of various standoff 
distances from 1 to 10 mm. 

Computed and observed main arc reverse polarity electric 
voltages in terms of various standoff distances from 1 
to 10 mm. 

Computed and observed crown widths in terms of various 
standoff distances from 1 ato 10 mm. 

Computed and observed root widths in terms of various 
standoff distances from 1 to 10 mm. 

Computed and observed average bead widths in terms of 
various standoff distances from 1 to 10 mm. 
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